Ordering dynamics of blue phases entails kinetic stabilization of amorphous networks 
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The cubic blue phases of Uquid crystals are fascinating and technologically promising examples of hierar- 
chically structured soft materials, comprising ordered networks of defect lines (disclinations) within a liquid 
crystalline matrix. We present the first large-scale simulations of their domain growth, starting from a blue 
phase nucleus within a supercooled isotropic or cholesteric background. The nucleated phase is thermody- 
namically stable; one expects its slow orderly growth, creating a bulk cubic. Instead, we find that the strong 
propensity to form disclinations drives the rapid disorderly growth of a metastable amorphous defect network. 
During this process the original nucleus is destroyed; re-emergence of the stable phase may therefore require 
a second nucleation step. Our findings suggest that blue phases exhibit hierarchical behavior in their ordering 
dynamics, to match that in their structure. 

PACS numbers: 61.30.Mp,64.60.qe 



The blue phases, BPI and BPII, of chiral nematic liquid 
crystals can each be viewed as an ordered network of topologi- 
cal defects (disclination lines), embedded within a liquid crys- 
talline matrix (a cholesteric) whose local molecular alignment 
axis rotates helically on moving through the sample. Such 
phases were long viewed as purely of academic interest, due 
to their extremely narrow window of thermodynamic stability 
(of order IK) LZI- This has changed completely with the re- 
cent development of new compounds showing stable BPs over 
a 50K interval, with fast switching between different states 
|2, 3 1. Thus BPs now offer a promising device technology not 
only for displays (H but also for laser applications lO [6l . To 
fully realize this potential requires an understanding of how 
chiral nematic materials switch between different structures. 
However, theoretical work on BPs has advanced relatively 
modestly since the late 1980s |7|, and so far there is almost 
no understanding of their phase-change kinetics. This is partly 
because of computational challenges which, despite pioneer- 
ing progress using small systems |[8l[9||9l, have prevented the 
simulation of supra-unit cell behavior in the ordered BPI and 
BPII, and ruled out realistic simulation of a third blue phase, 
the apparently amorphous | 7 1 BPIII. 

With the aid of supercomputers and a hybrid lattice Boltz- 
mann algorithm 1 1 1 we have recently overcome these difficul- 
ties, enabling us to address systems hundreds of times larger 
than the unit cell volume. This has allowed us to address 
elsewhere the motion of planar interfaces between compet- 
ing phases |2|. By the same methods, we address below for 
the first time the domain growth of the ordered BPs from an 
isolated nucleus of the stable phase. We find that the growth 
process is unexpectedly interrupted by the proliferation of a 
metastable BPIII-like structure. Thus the system seemingly 
opts for speed rather than efficiency, lowering its free energy 
rapidly by amorphous defect proliferation even though this 
creates barriers that prevent attainment of the global free en- 
ergy minimum - which could have been reached directly by 
a slower, orderly growth of the initial nucleus. This finding, 
which may have wide implications for phase kinetics in other 
multi-scale soft materials (see e.g. 1 12, 13 1), is sharply differ- 
ent from some other instances where metastable phases inter- 
vene in phase ordering ('Ostwald's rule') as we discuss below. 



Experimentally BPs are found adjacent to the cholesteric 
phase which has local nematic orientational order, along a di- 
rector that slowly rotates in space; in the cholesteric, there 
is a single rotation axis and the wavelength in this direction 
is called the pitch. (All these phases lie below the isotropic 
phase in temperature; see Fig.l for a phase diagram.) Un- 
like the cholesteric, the BPs have director modulation in all 
directions and comprise an interpenetrating lattice (BPI/II) or 
disordered network (BPIII) of so-called 'double twist cylin- 
ders' 121. These accommodate a higher degree of director 
twist than the simple cholesteric, at the price of creating a 
complementary network of topological defects (disclination 
lines) on which the nematic ordering vanishes; the resulting 
defect-lattice unit cells in BP I/II are shown in Fig.l. The 
blue phases are more stable than the simple cholesteric at high 
enough molecular chirality. 

It is known that, on the scale set by the isotropic-cholesteric 
transition, the free energy differences among BPs are rela- 
tively small. Nonetheless, per unit cell, these differences typ- 
ically remain large compared to the thermal energy, knT fT\. 
This suggests only a limited role for thermal noise in BPs (but 
see 1 10] for a contrary view). Thus one might expect strong 
hysteresis and metastability of phases, but this has not been 
reported experimentally |2l|3]lll[I31. It is possible |16 | that 
small nuclei of BPs persist on heating far into the isotropic 
phase, allowing hysteresis-free growth on reducing the tem- 
perature again. Whatever the source of the required nuclei, 
simulation studies of their subsequent growth can illuminate 
fundamental issues in phase transition dynamics, as exempli- 
fied by previous work on hard-sphere colloids I ITHTsI . 

Our hybrid lattice Boltzmann technique is summarized in 
the methods section and in supporting information (SI; Text 
SI), and detailed elsewhere LinSJ. Briefly, it marries a finite 
difference code for the advective relaxation of the order pa- 
rameter tensor Q(r) (with r spatial position) IH, governed by 
a suitable free energy functional J-'[Q] 13, to an efficient par- 
allel lattice Boltzmann (LB) code for a forced Navier- Stokes 
equation for momentum transport, like that used previously to 
simulate multiscale colloidal materials |6|. 

We have found in previous simulations 121 that a semi- 
infinite slab of stable BPI or BPII will invade an adjacent 
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region of supercooled isotropic or cholesteric phase. In the 
present work we extend these studies to the more reaHstic case 
of a locaHzed nucleus, comprising a few unit cells (typically 8) 
of the stable BP. This nucleus is prepared by excising it numer- 
ically from a periodic BPI or II, placing it within the chosen 
environment, and then allowing it to anneal (see Methods sec- 
tion). Annealing is done very close to the coexistence bound- 
ary between phases, where there is no tendency for the nucleus 
to grow; during annealing its surface reconstructs to achieve 
local equilibrium with the surrounding phase. We then quench 
the system, altering the reduced temperature r and/or dimen- 
sionless chirality k, (defined in supporting Text SI) so that the 
nucleated BP is now the equilibrium structure. In most of our 
simulations there is no thermal noise; adding modest amounts 
of noise makes quantitative but not qualitative differences to 
our findings (with one exception detailed below). 

Fig.2(a-c) shows the time evolution of a BPII seed in an 
isotropic environment. Depicted is an 'isosurface' at a fixed 
small magnitude of nematic order (defined through the largest 
eigenvalue < q < 1 of the order parameter tensor Q; see 
methods section). In a bulk BP this representation allows 
unambiguous imaging of the disclinations, each of which is 
wrapped by a thin tube of the isosurface. Here, since there 
is no nematic order in the external phase, this isosurface en- 
velops the growing nucleus; to expose the structure within, we 
omit in Fig. 2b the front section of the growing droplet. Al- 
though there is residual cubic anisotropy (smaller with ther- 
mal noise than without) the phase formed does not have the 
long range order of BPII; it is instead an amorphous, aperi- 
odic network. This disorder is also visible in the structure 
factor C(k) |[T9ll shown in Fig. 2d. By the end of the run, no 
trace remains of the initial nucleus. 

Fig. 3 shows a similar sequences for a BPI seed in a 
cholesteric matrix. The pattern of evolution is similar al- 
though the residual anisotropy is no longer cubic (Fig. 3d), 
as befits the lowered initial symmetry. The dynamics of the 
defect proliferation can be clarified further by using instead 
an extended rod-like nucleus of square cross section, whose 
long axis spans the periodic simulation box. When viewed 
along that axis (Fig.4) the noise-free growth of the disclina- 
tion pattern, though aperiodic, retains the symmetry of the 
initial state. Addition of thermal noise allows this symmetry 
to break and increases the rate at which order is lost (see sup- 
porting Fig.Sl); this is the expected consequence of adding 
weak noise to a time evolution governed by deterministic, but 
probably chaotic, dynamics. Finally, in Fig. 5 we show growth 
of a BPI seed in an isotropic environment. In this case, the 
quench parameters are the same as in Fig. 2, so that the nu- 
cleated phase BPI is metastable, with BPII stable. As before, 
an amorphous network forms in preference to either ordered 
phase. Supporting Videos S1-S4 show the complete time evo- 
lutions corresponding to Figs. 2-5. 

These results demonstrate that seeding a supercooled 
isotropic or cholesteric phase with a small BPI or BPII nucleus 
leads to growth, not of the stable phase, but of metastable, 
amorphous networks. Although clearly BPIII-like, these 
states probably do not belong to a metastable branch of BPIII 
itself (see supporting Text SI). A more important question. 



especially given the observed reversibility of the experimental 
phase diagram, is whether such amorphs can thereafter be- 
come ordered. We find that, at least for low chirality (k, <2), 
this transformation is not spontaneous: without noise, the sys- 
tem relaxes almost to a standstill, suggesting that a second 
nucleation is needed to reach equilibrium. We thus tried di- 
rectly re-inserting a BPI/II nucleus similar to the initial one; 
this disappears, indicating that the critical nucleus for this step 
is relatively large. At higher chirality, we see instead contin- 
uous ripening of the amorph towards of a long-range ordered 
phase (supporting Fig.S2). However, this is not BPI or II, but 
another phase (O^) long predicted by theory |7 | but not yet 
seen experimentally. More generally, our findings admit a va- 
riety of kinetic pathways that might connect the metastable 
amorph formed initially to BPI/II, although in some cases this 
final equilibration might not be observable on any experimen- 
tal timescale. 

At first sight, our results are reminiscent of Ostwald's rule 
of stages 1241 . whereby nucleation is asserted to proceed 
through a series of metastable phases, starting with the one 
whose free energy lies closest to (but below) that of the ini- 
tial bulk. There is however no simple link: Ostwald's rule 
is normally attributed to having a higher barrier for nucleat- 
ing the final structure directly |25|, whereas our simulations 
begin with a suitable nucleus already present. Moreover we 
have performed simulations with r and/or adjusted (sup- 
porting Tables SI, S2) to ensure that the initially nucleated or- 
dered BP now lies closest in free energy to the initial isotropic 
phase. An amorphous network still forms, now with free en- 
ergy lower than the nucleated phase, in contradiction to Ost- 
wald's rule. 

In conclusion, we have shown that the growth of a nucleus 
of an ordered BP within an isotropic or cholesteric matrix is 
generally interrupted by proliferation of a metastable amorph, 
despite the fact that this creates a barrier to formation of the 
equilibrium phase which could, however, have been created 
directly from the initial nucleus, whose structure it shares. 
Our results illustrate the subtle character of domain growth in 
soft materials exhibiting hierarchical ordering in equilibrium. 
For BPs, the organizational levels are local cholesteric order- 
ing; defect formation within the cholesteric matrix; and the 
long-range organization of the defects into a periodic lattice. 
Our results suggest a corresponding multi-stage time evolu- 
tion, which can easily get stuck at an intermediate level in the 
hierarchy. Once the amorphous defect network has formed, 
the energetics is already optimized at the first two levels; the 
driving force for creating an ordered BPI/II superstructure is 
then relatively weak, and the barriers high. Conversely, the fi- 
nal thermodynamic preference for that ordered structure, even 
if a nucleus of it is initially present, is apparently swept aside 
in the rush to create disclinations during the earlier defect- 
proliferation stage. These findings may have wide implica- 
tions for the phase ordering dynamics of multiscale soft mate- 
rials — increasingly relevant to nanotechnology, photonics, 
and display applications — of which the blue phases offer 
striking and elegant examples. 
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FIG. 1: Left and centre: disclination network of the crystalline Blue Phases BPI (left) and BPII (centre). The pictures show 8 unit cells, 2 
along every coordinate direction. Right: computed phase diagram in the plane spanned by chirality (k) and effective temperature (r). At high 
chiralities BPII becomes unstable relative to either the 05 structure (in theory) or the BPIII structure (in experiment); these phases are not 
shown here. The dotted arrows show the various quenches used in Figs. 2-5, from states close to a phase boundary to the final parameters at 
which domain growth occurs. 




FIG. 2: Domain growth of a BPII droplet in isotropic environment. The droplet was equilibrated near the BPII-isotropic boundary and then 
quenched to r = 0, At = 2, where BPII is the equilibrium phase. The pictures show isosurfaces (q = 0.13) of the scalar order parameter 
during equilibration, and at intermediate and late times (t = 4 x 10^; t = 7.6 x 10^; t = 6 x 10^) (a-c). (For typical materials the simulation 
timestep corresponds to ^ 1 ns.) Frame (d) gives a cut through the structure factor of the final state at A;^ = for wave vectors kx and ky in 
[—7v/2i, 7v/2i]. Here i is the discretization length of the LB lattice (about lOnm for typical materials). 



I. EQUATIONS OF MOTION 

The free energy functional |5 1, and the equations of motion 
m deriving from it, are presented in supplementary Text S 1 . 
As detailed there, the control parameters of the free energy 
functional are conventionally chosen as r, a reduced temper- 
ature, and a dimensionless chirality. We also report there 
our estimates of the free energy differences among BPs and 
specify the way in which thermal noise is introduced. 

II. QUENCH DETAILS 

For each quench shown in Figs. 2- 5, we report in supple- 
mentary Table S2 the free energy densities for the stable 
BPI/II phases and for the metastable BPIII-like end state. Ev- 
ery simulation comprised two separate parts, namely the cre- 
ation and equilibration of the initial configuration, followed 
by a quench in temperature and/or chirality. In order to gen- 
erate an equilibrated nucleus, we started from the analytically 



known structures | 7 1 that describe the infinite chirality limits 
of BPI and BPII and first evolved these for 2000 time steps at 
the chosen initial tz and r. This creates an equilibrium bulk 
phase at those parameters. Then we replaced the tensor or- 
der parameter at all sites outside of a cube comprising 8 BP 
unit cells, situated at the centre of the simulation box, with 
either a cholesteric or an isotropic configuration. This cre- 
ates a nucleus of BP crystal within the required phase, but 
the nucleus has sharply-cut surfaces that could have unrealis- 
tically high surface energy, perturbing the subsequent dynam- 
ics. This state was therefore equilibrated in a second stage for 
another 18000 time steps, maintaining r and k, at the chosen 
initial values, which lie just within the BP equilibrium region 
but very close to the coexistence line with the surrounding 
phase. By this process we allowed the surface structure and 
shape of the BP nucleus to equilibrate locally, without initi- 
ating significant domain growth. After the equilibration was 
finished, we performed a sudden quench to parameter values 
further away from the phase boundary, thereby initiating the 
domain growth. Note that for technical reasons (the require- 



4 




FIG. 3: Domain growth of a BPI droplet in cholestric environment. The droplet was equilibrated near the cholesteric-BPI boundary and then 
quenched tor = —0.5, At = 1.2, where BP I is the equilibrium phase. The pictures show isosurfaces (g^ = 0.2) of the scalar order parameter 
immediately after the quench at t = 2.4 x 10^ timesteps (a); and at two later stages, t = 5 x 10^ and t = 6 x 10^ (b and c). The top left 
half in (d) is a cut through the structure factor C(k) of the final state (c) along ky = for kx < kz in the interval [—n/2i,7v/2i], whereas the 
bottom right half depicts a cut along at /c^ = for kx > ky in the same range. 




FIG. 4: Domain growth of a BPI rod in cholestric environment: The nucleus was equilibrated close to the BPI-isotropic boundary and then 
quenched to r = — 0.5, A€ = 1.35. The pictures show isosurfaces (q = 0.2) of the scalar order parameter immediately after the quench at 
t = 2.4 X lO'^ timesteps (a), and at two later stages t = 2 x 10^ and t = 3 x 10^ (b and c). (d): The structure factor of state (c) at /c^ = for 
wavevectors kx (horizontal axis) and kz (vertical axis) in the interval [—7T/2i,7T/2£]. 



ment of fitting a whole number of unit cells into the simulated 
domain), when nucleation is done within a cholesteric matrix, 
the cholesteric pitch, post-quench, is not the optimal one. The 
same would apply in generic laboratory conditions involving 
samples of fixed geometry where a pure temperature quench 
will typically alter qq as well as both r and k,. 



lattice of 128^ sites capable of containing around 500 BPI/II 
unit cells. For typical cholesteric materials, the LB lattice 
spacing corresponds to about 10 nm, and the timestep about 1 
ns (see supplementary Text SI). Similar parameter mappings 
arise in LB simulations of hierarchical colloidal materials (61. 



III. STRUCTURAL VISUALIZATION 



The openDX software package was applied to the simulated 
q{T) data to create the constant-g isosurfaces shown in Figs.l- 
5. Choosing a small finite q renders the disclination lines 
as tubes without topological ambiguity; the visually optimal 
choice of q depends on the structure being viewed. The struc- 
ture factor C(k) is defined as |^(k)p with g'(k) the Fourier- 
transform of q{r) computed using routines from the FFTW 
library. 

IV. COMPUTATIONAL RESOURCES 

A typical computational run time was 5-10 hours on 512 
quad-core nodes of an IBM Blue Gene/P system, using a LB 



Acknowledgments 



We acknowledge EPSRC (Grants EP/ E045316, EP/ 
E030173, EP/ F054750 and EP/ C536452) for funding and 
supercomputer time. We also acknowledge a DOE INCITE 
2009 grant entitled 'Large scale condensed matter and fluid 
dynamics simulations' at Argonne Leadership Computing Fa- 
cility, supported by the US Department of Energy under con- 
tract DE-AC02-06CH11357. MEC holds a Royal Society Re- 
search Professorship. 



5 




FIG. 5: Domain growth of a BPI droplet in an isotropic environment. The droplet was equilibtated close to the isotropic-BPI boundary and 
then quenched to r = 0, k = 2, where BPII is stable. (See phase diagram in Fig.l.) The pictures show isosurfaces (q = 0.13) of the scalar 
order parameter: during the equilibratio and at intermediate and late times (at t = 4 x 10^ ;t = 1.2 x 10^; t = 6 x 10^) (a-c). In picture (b) 
the growing domain has been cut at ^ = 32 to reveal the internal structure. A cut at A:^ = through the stucture factor C(k) of the final state 
is depicted in (d). The section shown corresponds to wavevectors kx (horizontal axis) and ky (vertical axis) in the interval [—7v/2i,7v/2i]. 
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Appendix A: Supporting Text 

Here we provide further information on the free energy 
functional and the control parameters r, k, (section 1); on the 
equations of motion (section 2); on the conversion of simu- 
lation parameters to physical units (section 3); on the intro- 
duction of thermal noise (section 4); on the parameter values 
for our quenches (section 5); and on the end- state free ener- 
gies in comparison with various stable and metastable phases 
(section 6). 

1. Free energy functional, chirality and reduced temperature 

The thermodynamics of cholesteric blue phases can be de- 
scribed via a Landau-de Gennes free energy functional 
which in turn is an integral over space of a free energy density 



^[Q] 



d^r/(Q(r)). 



(Al) 



The free energy density / = /(Q) may be expanded in pow- 
ers of the order parameter Q and its gradients; Q is a trace- 
less and symmetric tensor. The largest eigenvalue q and cor- 
responding eigendirection (n) respectively describe the local 
strength and major orientation axis of molecular order. The Q 
tensor theory, rather than a theory based solely on the direc- 
tor field n(r), allows treatment of disclinations (defect lines) 
in whose cores n is undefined; in blue phases, disclinations 
organise into regular or amorphous networks. 
The free energy density we use, following Q, is: 



/(Q) = ^(i-|)« 



2 



4 



K 



K 



Here repeated indices are summed over and stands for 
Qa(3Qa(3 ^tc. The first three terms are a bulk free energy 
density whose overall scale is set by Aq (discussed further be- 
low); 7 is a control parameter, related to reduced temperature. 
Varying the latter in the absence of chiral terms (qo = 0) gives 
an isotropic-nematic transition at 7 = 2.7 with a mean-field 
spinodal instability at 7 = 3. 

The rest of the free energy in Eq |A2| describes distortions 
of the order parameter field. As is conventional (5] [3 we as- 
sume that splay, bend and twist deformations of the director 
are equally costly; K is then the one elastic constant that re- 
mains. The parameter qo is related via qo = 27r/po to the 
pitch length, po, describing one full turn of the director in the 
cholesteric phase. In BPs one observes that the pitch length 
(still defined locally by the spatial rotation rate of n) slightly 
increases on entering the BP from the cholesteric phase. To 
account for this, a 'redshift' factor r is introduced |9 1 whose 
variation effectively allows free adjustment of the BP lattice 



parameter, A ^ A/r. To avoid changing the size of the simu- 
lation box, redshifting is performed in practice by an equiv- 
alent rescaling of the pitch parameter and elastic constant, 
^0 Qo/'^ and K K r'^ . In simulations aimed solely at 
free energy minimization, it is legitimate to make r a dynamic 
parameter and update it on the fly to achieve this 1 9 1 . We do 
not allow this during our domain growth runs, but do relax r at 
the end of selected runs as part of our free energy comparison 
(see Section 6). 

The phase diagram of blue phases (if thermal fluctuations 
can be neglected ifTOl ) depends on just two dimensionless 
numbers, which are commonly referred to as the chiral- 
ity, and r, the reduced temperature |7|. In terms of the above 
parameters, these are: 



27(1 - 7/3) 



Aoi 



(A3) 
(A4) 



If the free energy density Eq |A2| is made dimensionless, r ap- 
pears as pref actor of the term quadratic in Q, whereas n quan- 
tifies the ratio between bulk and gradient free energy terms. 



2. Equations of motion 

A framework for the dynamics of liquid crystals is provided 
by the Beris-Edwards model |4|, in which the time evolution 
of the tensor order parameter obeys 

{dt + uo^do^) Q - S(W, Q) = TH + C (A5) 
where C is a noise term discussed in section 4. In the absence 



of flow, Eq |A5 1 describes relaxation towards equilibrium, with 
a rotational diffusion constant F, driven by the molecular field 
H. The latter is given by | 4 1 



I [ST 
SQ 3 ^ \SQ 



(A6) 



The tensor S in Eq |A5 [ couples the order parameter to the sym- 
metric and antisymmetric parts of the velocity gradient tensor 

Waf3 = df3Ua, defined as 



A = (W + W^)/2; n = (W -W^)/2. 
This coupling term reads explicitly 



(A7) 



S(W,Q) = (^A + l^)(Q+^) + (Q + ^)(^A-l^) 



2^(Q + ^)Tr(QW). 



(A8) 



Here ^ is a material-dependent 'tumbling parameter' that con- 
trols the relative importance of rotational and elongational 
flow for molecular alignment. We choose ^ = 0.7, within the 
'flow aligning' regime for which molecules align at a fixed 
angle (the Leslie angle) to the flow direction in weak simple 
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shear O. (For 'flow tumbling' materials, the director instead 
rotates continuously p^.^Sll.) 

The momentum evolution obeys a Navier-Stokes equation 
driven by the divergence of a stress tensor T,: 



(A9) 



This pressure tensor is in general asymmetric and includes 
both viscous and thermodynamic components 



(AlO) 



Note that within the lattice Boltzmann flow solver, the 
isotropic pressure and viscous terms are managed directly by 
the solver (as in a simple Newtonian fluid, of viscosity r]) 
whereas the divergence of the remaining terms is treated as 
a local force density on that fluid. 



3. Parameter mapping to physical units 

Here we describe how simulation parameters are related to 
physical quantities in real BP materials. In order to get from 
simulation to physical units, we need to calibrate scales of 
mass, length, and time (or equivalently length, energy and 
time). We follow a methodology similar to that of | 8 1. 

First we define a set of LB units (LBU) in which the lattice 
parameter i, the time step At, and a reference fluid mass den- 
sity po are all set to unity. This is the set of units in which our 
algorithm is actually written. The first two of these param- 
eters directly connect to observables in the simulations and 
we show below how to map these onto experiments. The fluid 
density, however, enters differently. So long as fluid inertia re- 
mains negligible (low Reynolds number. Re = pVA/rj, where 

A are typical length scales and velocities of any flow) the 
physics observed will not depend on the actual mass density p 
of either the physical or the simulated system. Since LB uses 
inertia to update the fluid velocities, it improves efficiency to 
use a density that causes Re to be several orders of magnitude 
larger than in experiments; so long as Re remains small (say. 
Re < 0.1), no harm is done |11|. Such parameter steering is 
helped by allowing p ^ po within the code. 

We now turn to the calibrations of length, energy and time. 
The length scale calibration is straightforward, and fixed by 
the cholesteric pitch po, which is typically in the 100-500 nm 
range | 7 |. More precisely, in our simulations we set the unit 
cell of BPI/II to be 16 LBU; this gives good resolution without 
wasting resource. Therefore the LBU length unit (one lattice 
site) corresponds to, say, lOnm in physical space. 

To get an energy scale, we use the measurements cited in 
Appendix D of | 7 |, which suggest 



Here a, 6 are parameters defined in fTl, and the required ratio 
is expressed in terms of our chosen parameters (whose com- 
parison with |7| shows that b^/6 = and a = Aqj/A). 
From this relation, given that ^ c::::^ 3, we obtain that Aq — \0 
Pa. On the other hand, our simulations use a value of Aq ~ 
0.01 LBU. This requires that the LB unit of stress is about 
10^ Pa in SI units. We next use this to calibrate the LB time 
unit and then crosscheck that the resulting fluid density gives 
acceptable Reynolds numbers. 

For the timescale calibration we use the following formula 
(HIHl which relates the so-called 'rotational viscosity' 71 (de- 
fined by the relaxation equation jidn/dt = —SF/Sn) for 
the director field in a well- aligned nematic, to the ordering 
strength q and the order parameter mobility F: 



71 



r 



(A12) 



a 



27 

2A^ 



5 X 10"^J"^m^ 



(All) 



In our simulations, we choose F = 0.3 and also select the 
thermodynamic parameters to give q 0.3. Therefore 71 :^ 
0.6 LBU. For real materials, 71 lies usually in range 10~^ — 
10~^ Pa s 1 5 1; we choose for definiteness 71 = 0.06 Pa s 
= 0.6 LBU. Given the previous result for stress, this requires 
that the LB time unit (one algorithmic timestep) equates to 
10-^ s. 

Note that in our simulations we use a fluid viscosity r] 
1 — 2 LBU = 0.1 — 0.2 Pa s (using the parameter mapping just 
established). This is sensible (if somewhat low) for a molecu- 
lar nematogen in the isotropic phase. (The effective viscosity 
in ordered phases is of course higher, but the coupling to the 
order parameter handles this.) Similarly, we adopt elastic con- 
stant values K 0.01 LBU = 10~^^ N, corresponding to a 
Frank elastic constant K/2 of 50 pN which is again sensible 
15]. 

Finally, we need to cross-check the fluid density. We have 
1 kg m-3 = 1 Pa s^ m-2 = iQ-s+is-ie lbU = 10"^ 
LBU. Thus the reference density po equates to a fluid density 
p = 10^ kg m~^, roughly a thousand times larger than exper- 
imental values. As explained previously, however, this makes 
no difference so long as the Reynolds number Re = pVA / r] 
remains small enough. This dimensionless number can be 
evaluated directly in lattice units. We have p and r] of order 
unity, and set A ~ 16 (a BP unit cell). Observing that typical 
velocities arising in our simulations are around 10 ~^ LBU we 
get Re < 10~^. Even allowing for higher peak velocities and 
larger A in some materials, this is safely small ifTTIl . 

To summarize the above, our simulations faithfully repre- 
sent experimentally realisable BP-forming materials, subject 
to the interpretation of the LB units for length, time, and en- 
ergy density are close to 10 nm, 1 ns, and 100 MPa respec- 
tively. 



4. Adding noise 

With these choices, the typical difference in free energy 
density between say BPI and BPII is 10"^ - 10"^ LBU (see 
Table S2) or about 100 - 1000 Pa. Even for a rather small BP 
lattice constant of 160 nm, this corresponds to a free energy 
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difference per unit cell of 100 — 1000 ksT, much larger than 
any thermal energy. Such unit-cell level differences represent 
a reasonable estimate of the barrier to topological reconnec- 
tions within an evolving BP structure, and appear to preclude 
any significant role for thermal noise in the domain growth 
process. Nonetheless it has been argued theoretically ifTOl that 
entropy plays an important role in BP thermodynamics and 
presumably also therefore, domain growth. (This certainly 
becomes more likely for atypically small Aq and/or small unit 
cells.) We therefore repeated selected simulations with ther- 
mal noise present. 

We choose to add noise to the order parameter sector only, 
via the term in Eq |A5[ In principle, noise can also be 
added to the fluid mechanical sector |12]; this will create 
an additional conserved diffusion of the order parameter (a 
kind of Brownian motion) which should however be negligi- 
ble, at least at large length scales, compared to the local non- 
conservative relaxation embodied in Eq |A5[ The fluctuation- 
dissipation theorem then fixes the variance of the noise in 
Eq|X5]as 

{Ca^{r,t)C^s{r,t')) = 2kBTrPa^^sS{r - r')S{t - t'). 

(A13) 

Here P is a tensor that projects the order parameter into inde- 
pendent components (respecting its symmetry and traceless- 
ness) and allows only diagonal correlations in that space. 

By using the parameter mapping detailed above, we can 
find the value of ksT = 4 x 10~^^Jin simulation units as 
4 X 10~^ LBU. However, to allow for variation in material 
parameters, such as the BP unit cell size, we have studied a 
range of noise temperatures in the range 10~^-10~^ LBU. 
At the upper end of this range we can observe vibrant ther- 
mal fluctuations of the disclination network, associated with 
some shifting of phase boundaries. This merits further study 
in view of the claims of 1 10] that fluctuations can signigicantly 
alter the equilibrium physics of blue phases, although a full 
exploration of such effects lie beyond the scope of the cur- 
rent work. However, in terms of domain growth dynamics, 
we have found the main effects of adding modest amounts of 
thermal noise to be quantitative rather than qualitative. (Here 
'modest' means /c^T < 10~^ LBU, as might be relevant to 
BPs with lattice parameters > 300 nm rather than the 160nm 
value used in the parameter mapping in section 3.) An excep- 
tion was already discussed in relation to Fig.Sl, which shows 
that quite low noise levels can disrupt the orderly but aperiodic 
growth that would otherwise arise from a rod-like nucleus. 



5. Quench details 

Quenches were performed as described in the Methods sec- 
tion of the main text. Supporting Table S 1 specifies the com- 
plete simulation parameters (all in LBU) for each of the runs 
presented in Figures 2-5 of the main text and Figures S1-S2. 



6. Free energy comparisons 

Supporting Table S2 compares the free energy densities 
of three crystalline blue phases (BPI, BPII, 05) with the 
cholesteric phase and the end- state amorphous networks for 
the final r, k, values appropriate to all runs reported in Figs. 2-5 
and Figs.Sl,S2. (The noise-free isotropic phase has free en- 
ergy density zero in all cases.) Note that 05, as long predicted 
theoretically Q, becomes the most stable ordered phase at 
high chirality. However this phase is not seen experimentally 
and its relative stability might be a consequence of the one 
elastic constant approximation or some other shortcoming of 
the free energy, Eq |A2[ 

For end-state amorphs nucleated within an isotropic phase, 
two free energy values are reported. The higher one is the 
direct result of the quench, performed at fixed redshift (as is 
appropriate when simulating directly the equations of motion 
presented above). The lower value is for a annealing protocol 
whereby the redshift is released at the end of the main part of 
the simulation, either when the disclination network has filled 
the simulation box and rearrangements have come to a virtual 
standstill (this was done for r = 0, = 2) or when the discli- 
nation network first collides with its own periodic images (this 
was done for r = 0, = 3). This 'devil's advocate' anneal- 
ing schedule finds the best possible free energy among states 
with topology close to the amorph, regardless of whether such 
states were dynamically accessible from the initial nucleus. 

For < 2, the final redshift release makes no qualitative 
difference to the observed topology; moreover the free en- 
ergies found after this procedure for (k, = 2) still lie above 
those of the ordered 'target' structure, as they do for all runs 
with K <2 done in the absence of a final redshift-release step. 
This confirms that the amorphous end- state is metastable, as 
claimed in the main text. 

The quenches at r = 0, = 3 offer a test of Ostwald's rule 
of stages; for these parameters, the free energy of the final 
amorphous network (with or without redshift release, the lat- 
ter in bold) lie further in free energy from the initial isotropic 
phase (/ = 0) than the structure that was nucleated. When 
this is BPI, so do all remaining ordered phases. Thus the nu- 
cleus is now made of the phase whose free energy lies closest 
below that of the initial bulk; by Ostwald's rule, this should be 
the first phase to grow. Despite this, an amorphous network is 
formed, and the initial nucleus disappears - just as it did when 
the nucleated phase was the stable one. Therefore our results 
have no explanation in terms of Ostwald's rule; indeed, for 
this system, they disprove it. 

In these high chirality runs, the final release of the redshift 
noticably changes the dynamics, accelerating the continuous 
ripening of the amorphous network towards the equilibrium 
05 phase. No redshift release was performed for runs in- 
volving nuclei within a cholesteric matrix (r = — 0.5, = 
1.2, 1.35) due to the much longer simulation times required 
in this case. Indeed, for runs (F2,F3) the defect network did 
not fill the box by the end of the run and the quoted free en- 
ergy densities in these cases are upper bounds. However, cor- 
recting for such volume-fraction effects shifts the values for 
the amorphs down by less than one percent, confirming their 
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metastability. 

The amorphous networks we report in this work appear 
closely related to BPIII, an equilibrium blue phase, believed 
to be amorphous, which is found experimentally at high chi- 
ralities. We address elsewhere (manuscript in preparation) the 
question of whether the chosen free energy density, Eq |A2[ 
can indeed predict a stable BPIII phase in that regime. Of 
significance here is the fact that it does predict a metastable 
BPIII phase for r = and n = 2^?> (relevant to Figs. 4, 5, 
and S2). To establish this, we have created candidate BPIII 
structures by evolving an initial state consisting of randomly 
oriented and positioned double twist cylinders (DTCs), em- 
bedded in a cholesteric matrix. On suitable annealing (with 
redshift enabled), these relax to form amorphous networks. 



metastable relative to BPI/II, that on visual inspection look 
very similar to our end- state structures. However, for r = 
0, = 2, 3, the DTC-based networks have a lower free energy 
(/ = lOV = -3.104, -1.144) than the amorphs formed in 
our nucleation runs with BPI/II. (Except for the case of the 
BPII-nucleated amorph dX n = 3, this remains true even if 
we violate the true dynamics by releasing the redshift.) Thus 
our end- state amorphs are not the disordered network phase 
of lowest free energy in any of these cases. Arguably however 
BPIII is that phase: certainly it must be so at higher chirality, 
where it is thermodynamically stable. Thus we do not desig- 
nate our amorphs as belonging directly to a metastable branch 
of BPIII. 
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TABLE I: [Table SI] Simulation parameters in lattice units: temperature Teg, chirality /^eg, bulk free energy constant Ao,eg and 7eg at the 
equilibration point, temperature r, chirality bulk free energy constant Aq and 7 after the quench, redshift r as used during the equilibration 
and the main part of the simulation, elastic constant K, and noise temperature ksT. For all runs the pitch qo — 7r/16; rotational diffusion 
constant F = 0.3; the tumbling parameter ^ = 0.7; the viscosity 77 = 5/3, and the fluid density p = 2. 
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TABLE II: [Table S2] Free energy densities /(r) = 10^ f{r) in LBU for the ordered blue phases (BPI, BPII and 05) and the cholesteric 
phase (CH) and for the amorphous networks grown from a BPI-droplet in isotropic (BPI-ISO-D) or cholesteric environment (BPI-CH-D), 
a BPII-droplet in isotropic environment (BPII-ISO-D) or from a BPI-rod in cholesteric environment with (BPI-CH-RN) and without noise 
(BPI-CH-R). The free energy density /(r) of each phase was attained by direct evolution of the initial state at fixed redshift r; f{r), r' are 
corresponding values after redshift release. Numbers in bold provide tests of Ostwald's rule(see supplementary Text SI, section 6); italics 
identify the equilibrium phase at given r, n. 
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FIG. 6: [Fig.Sl] Domain growth of a BPI rod in cholesteric environment including noise in the order parameter equation of motion. The 
noise strength was set to ksT = 8.33 • 10~^ LBU. The pictures show isosurfaces of the order parameter at t = 2.4 x 10^ timesteps (a); 
t = 1.2 X 10^ (b); t — 2 X 10^ (c); and t = 2.9 x 10^ (d). The equihbration and quench were performed in exactly the same way as in the 
simulation of the BPI rod without noise (Figure 4 in main text). Note that from about t = 10^ onwards, the disclination network ceases loses 
the symmetry of the initial state. 





FIG. 7: [Fig.S2] Domain growth of a BPI droplet in isotropic environment at high chirality. The droplet was equilibrated near the isotropic-BPI 
boundary at r = 0.95, k = 0.4 and then quenched to r = 0, At = 3. At this point 05 is the thermodynamically stable phase. The pictures 
show the isosurface {q — 0.12) of the scalar order parameter during the equilibration (a), at an intermediate stage (b) and in the final state of 
the simulation (c). In (b) the isosurface has been cut at 2/ = 32 to reveal the internal structure. In (d) a cut through the structure factor of the 
final state (c) along kz — for kx < ky on the interval [— 7r/2^, 7r/2^] is depicted. Although the Bragg-peak pattern appears rather blurred, 
the cubic signature of 05 is clearly recognizable. 



